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Abstract 

It is known, but perhaps not well-known, that when the mortality 
is assumed to be of Gompertz-Makeham-type, the expected remain- 
ing life-length and the commutation functions used for calculating the 
expected values of various types of life insurances can be expressed 
with an incomplete gamma function with a negative shape parame- 
ter. This is not of much use if ones software cannot calculate these 
values. The aim of this note is to show that one can express the 
commutation functions using only the exponential function, the (or- 
dinary) gamma function and the gamma distribution function, which 
are all implemented in common statistical and spreadsheet software. 
This eliminates the need to evaluate the commutation functions and 
expected remaining life-length with numerical integration. 

Keywords: Gompertz, Makeham, commutation functions, continuous 
compounding, analytical expression. 

We will assume that individual lifes are distributed according to the Gomp- 
ertz-Makeham distribution with distribution function F{x; a, (3,'~f) = 1 — 
l{x; a, (3, 7) and survival function l{x; a, (3, 7) = exp{—ax — ^{e^^ ~ ^)} cor- 
responding to the mortality rate /i(a;; a,/5,7) = —-^l{x; 7) = a + Pe^^ 
at age x. We assume a fixed continuously compounded interest rate 6. Let 
(ix = a,x{a, l3,'y,6) be the present value of a life-long annuity, continuously 
paid at rate 1, to a person at age x, and ea;(a,/3, 7) the expected remaining 
life- length for a person at age x. The main result is the following. 
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Theorem 1. 

eo{a + S, Pe'^^ where (1) 
i (l - {^r^'e^^^ni - ^) [1 - G(f ; 1 - ^, 1)] ) , z. e. (2) 

a+S \ \ 7 / 

xr(l-^)[l-G(^;l-^,l)]), (3) 

where T{ri) is the gamma function: T{ri) = y^~^e~^dy, and G{z;ri,l) is 
the distribution function of a gamma distributed random variable with shape 
parameter rj and scale parameter 1: G{z;ri, 1) = y^^^e'^dy. 

Remark 1. Even though equation ([3]) may look unwieldy, it only consists of 
functions that are readily available in common statistical and spreadsheet 
software, e.g. in R and Microsoft Excel. Thus, no numerical integration is 
necessary to evaluate these expressions. Note, however, the findings of Yalta 
(2008). 

Remark 2. The special case with zero mortality, i.e. a = (3 = (and 7 arbi- 
trary), yields a2,.(0, 0, 7, 5) = eo(5, 0,7) = 1/6 as expected for a continuously 
compounded perpetuity. With a > and /? = 0, the life-length is exponen- 
tially distributed and there is no ageing. Then a^^a, 0, 7, 6) = eo{a+6, 0, 7) = 
l/{a + 6). The general case has dx{a, /5, 7, 6) = -£^{1 — [•••]) and one can 
thus interpret the ellipsis part [. . . ] as the effect of ageing on the value of the 
annuity. 

With 6 = 0, the present value of the annuity equals the nominal value, and 
since the annuity is paid continuously at rate 1, this nominal value is equal 
to the remaining life-length. Thus, we arrive at the following result, which 
will also be reached through another route in the proof of Theorem 1. 

Corollary 1. e^.(a,/3,7) = a^(a,/3,7,0) = eo(a, /3e'^'^', 7). 

In order to calculate the present values of different life insurance contracts 
one usually defines the commutation functions 

D{x) = /(x)e-^^ 

POO 

N{x) = / D{y)dy, 

J X 

POO 

M{x) = / i^{y)D{y)dy. 



a^{a, (3,^,6) = 
eo(a,/5,7) = 

ax{a,P,j,6) = 
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The corresponding commutation functions with twice the interest rate are 
also useful when one wants to calculate not only the expected value of dif- 
ferent insurances, but also their variance, see Andersson (2005). In general 
M{x) = D{x) — 6N{x), and with Gompertz-Makeham mortality D{x) = 
l{x]a + 5, /9, 7), which is easy to evaluate numerically. The only "difficult" 
function is N{x). However, a standard result is that = N{x)/D{x), so by 
Theorem 1 we get a tractable expression also for N{x): 

Corollary 2. N{x) = D{x)eo{a + 6, f3e^^,'y), where D{x) = /(x; a + 5, /5, 7). 
It remains to prove Theorem 1. 

Proof of Theorem 1. The incomplete gamma function is defined as T{ri, z) = 
IT y^~^^~^^y- analogy with the gamma distribution we call 77 the shape 
parameter. The incomplete gamma function is closely related to the gamma 
distribution function: 

T{ri,z) 



G{z- ri, 1) = — - / y^-^e-^dy = 1 - — - / y^-^e-^dy = 1 - 
or, equivalently, 

Tir^,z) = Tir^){l-G{z;r^,l)), (4) 
for T) > 0. By partial integration, 



POO 1 1 POO 

r(r7, z) = / y^-'e-ydy = - [y\-y] ^ + - / y^^e-^dy 



- (T{ri + l,z)-z''e-') . 



Combining this with (jlj) yields 

r(r/, z) = (^"e-^ - r(r7 + 1) [1 - G{z- r/ + 1, 1)] ) , (5) 
which we will use shortly. 

Following Andersson (2005), we express Cq with an incomplete gamma func- 
tion: 

eo(a,/5, 7)= / l{t; a, P,'y)dt = / e~°'^~^^'^^ ~^^dt (6) 
Jo Jo 



= - / y---'e-ydy 

7 J V7/ 7 Ji3/j 



7/ 7 V 7 7 



3 



We continue with the expected remaining hfe-length at age x: 



Jo i'{x;a,p,'y) Jo 

and by comparing with we arrive at the conclusion of Corollary[T) e^ia, P, 7) 
eo(a,/?e'^'',7). 

By standard arguments, 



00 



Jo Ha;;a,/3,7) Jo 

and by once again comparing with equation ([T]) in the theorem is proved. 

Equation ([2]) is obtained by using ([5]) on (171) to rewrite the incomplete gamma 
function in the latter equation. Equation ([3]) follows immediately from ([T]) 
and ([2]). This concludes the proof of Theorem 1. □ 

Remark 3. Note that for typical values for the parameters, e.g. a = 0.001, j3 = 
0.000012,7 = 0.101314 and 6 = 0.026559 from Andersson (2005), we have in 
equation a shape parameter 1 — (a + 5)/7 = 0.727984 which is positive, 
just as shape parameters for gamma distributions should be. If one has to 
ones disposal software that can calculate incomplete gamma functions with 
negative shape parameters, then one can use equation ([7]) straight away, see 
e.g. Mingari Scarpello et al. (2006). Since the effect of ageing on the value of 
an annuity, cf. Remark [21 is not easily seen from ([71), the formulation in the 
theorem is nevertheless of some independent theoretical interest, regardless 
of whether it is easy to evaluate functions with negative shape parameters 
or not. 
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